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Abstract. We perform Monte-Carlo simulations to study the Bernoulli (p) bond 
percolation on the enhanced binary tree which belongs to the class of nonamenable 
graphs with one end. Our numerical results show that the system has two distinct 
percolation thresholds pci and pc2- The mean cluster size diverges as p approaching 
to Pel from below. The system is critical at all the points in the intermediate phase 
(pci < p < Pc2) and there exist infinitely many infinite clusters. In this phase the 
corresponding fractal exponent continuously increases with p from zero to unity. Above 
Pc2 the system has a unique infinite cluster. 
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Figure 1. (left) Schematic diagram of the enhanced binary tree with 6 generations. 
Tlie danghng bond on left and right side is connected in the periodic boundary 
condition, (right) Another presentation of EBT (thick-black line) and its dual lattice 
(thin-red line). For easiness to see, the root vertex is not displayed for both the EBT 
and the dual lattice. 

1. Introduction 

Geometry of a space, which is characterized by dimensionality and topology, is a very 
important factor for the collective dynamics and statics of the elements embedded on it. 
It takes special significance in critical phenomena where correlation length diverges. 
Critical phenomena have been one of the main subjects of statistical physics and 
intensively studied over various systems such as interacting spin models, map cellular- 
automata, percolation, etc. In recent years, more and more attention has been paid to 
critical phenomena on a class of graphs which have quite different properties with those 
on the well-studied Euclidean lattices. One representative example is the class, so-called 
blue(almost) transitive nonamenable graphs (NAGs) P,[2], e.g., Cayley trees, hyperbolic 
lattices and enhanced trees [3] (see figure [1] for the binary case). Here "transitive" 
graph means that all vertices on it play exactly the same role in the large size limit. 
The fundamental property of transitive NAGs is that the number of vertices in a ball 
centered at the origin increases exponentially with the radius. Recent studies have 
reported that this property crucially affects the collective dynamics to make it differ 
much from those on Euclidean lattices [Il[2lll,[5l[6l[7l[8l[9l[3l[l0]. Although it seems that 
such graphs are non-realistic and studies on them are only for theoretical interests, their 
distinctive features, so-called small-world propeitj and a non-zero clustering coefficient, 
are common with complex networks, which are widely found in the social and natural 
systems [TTl [T2]. For example, percolation problems on complex networks are related 
to error and attack tolerance and cascade failures of real world systems [13]. Instead 
of extensive studies the interplays between network structures and critical phenomena 
is still unclear. Transitive NAGs has large advantages in its simplicity, e.g., regularity 
of the lattice, compared with usually studied complex networks. Studies on them has 
have a possibility to give a breakthrough on this complex problem. 

The number of ends, which indicates a kind of vulnerability of graphs, is an 
important concept for NAGs. The number of ends e{G) of an infinite graph G is 
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defined as the supremum of the number of infinite connected components in G\S, 
where G\S is the graph obtained from G by removing arbitrary finite subgraph S and 
the bonds incident to those. While trees have infinitely many ends, hyperbohc lattices 
and enhanced trees have only one end. Recent analytical studies have indicated that 
the Bernoulli bond percolation on NAGs with one end exhibits a new type of phase 
transition, so-called multiple phase transition (MPT), which takes three distinct phases 
according to open bond probability p as follows [H [2] : 

• Nonpercolating phase: there is no infinite cluster for < p < Pci{G). 

• Intermediate phase: there are infinitely many infinite clusters for Pci{G) < p < 

Vc2{G). 

• Percolating phase: there is a unique infinite cluster for Pc2{G) < p < 1. 

Here infinite cluster means a cluster whose mass diverges with system size as A^'^ with 
< ip < 1. As well known for the Euclidean lattices, any amenable graph, e.g., can have 
at most only one infinite cluster, i.e., Pd = Pc2 By contrast there are always plural 
infinite clusters in the percolating phases of trees. This is a consequence of infinitely 
many ends, i.e., trees have no long cyclic paths and split easily. Both of an exponential 
growth and sufficient amount of loops are necessary for graphs to exhibit three distinct 
phases. Indeed, Benjamini and Schramm proved the existence of a MPT on any planar 
transitive NAG with one end [10]. In addition, it is proved that the transition at Pd 
belongs to the mean-field universality class fL5\ . 

Previous studies of MPTs have been restricted in the analytic way based on 
probability theory and we lack numerical research, which will bring us the quantitative 
indication of MPTs. Especially finite size effect is hardly investigated by the analytic 
approach. This effect is crrucial to argue the NAGs and the realistic network, where 
the typical length scale of the graphs often grows as a logarithmic function of the total 
number of vertices A^ and is still small for considerably large A^. In this paper, we 
investigate a MPT by Monte-Carlo simulations taking the enhanced binary tree (EBT) 
as an example. Our numerical results actually give two distinct critical probabilities pd 
and Pc2- AH the points in the intermediate phase {pd < p < Pc2) are critical and there 
exist infinitely many infinite clusters . The corresponding fractal exponent continuously 
increases with p from to 1 and bridges a gap between non-percolating and percolating 
phases. 

2. Fundermental mechanism of the two-stage transition 

In this section we provide a basic picture of MPTs on NAGs. Although we consider the 
EBT as an example, the concept is not specified to it. 

At first we define an enhanced binary tree with L generations. We prepare vertices 
Vn,m for n = 0, 1,---,L — 1 and m = 0, 1, ■ ■ ■ , 2" — 1. The total number of vertices 
A^ equals 2^ — 1, i.e., L ^ log2 A^ . Inter-generation (radial) bonds are supposed 
between Vn,m and Vn+i,2m and between Vn,m and Vn+i,2m+i- These bonds yield a binary 
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tree. Furthermore intra-generation (circumferential) bonds are added between Vn.m and 
Vn,m+ii which means "enhanced" . All of inter- and intra-generation bonds are open with 
homogeneous probability p to yield connected clusters. EBT is not exactly transitive 
but almost transitive since the root vertex vq^ is special and the center can be defined. 
But this presumably does not affect the essential property of the critical phenomena. 

The following scenario is based on the natural assumption that the connectedness 
(correlation) function; the probability that a vertex of the £-th generation is connected 
to the root clster, that fo,o belongs to, exponentially decays wiht £ as 

Co{^) oc 2-^/«. (1) 

Here ^ means a correlation length along the generation (radial) direction and it 
monotonically increases with bond opening probability p. The averaged mass of the 
root cluster (sq) is calculated by summing up the connectedness function as 

L-l I ^ 

(so) oc y 2^Co(£) = where a = 2^~^/^ (2) 

e=o 

and 2^ is the number of vertices at the i-th generation taking a role of surface area factor. 
For < 1 and a < 1, (sq) converges to finite value in the thermodynamic limit L — > oo. 
On the other hand (sq) diverges with L for ^ > 1 and a > 1. It is important that the 
mass of cluster diverges despite of finite correlation length owing to the exponentially 
diverging surface area. This is an exact mechanism for the percolation transition on the 
binary Cayley tree, where Co{i) = p^ and C, = l/log2(l/p). We can define the critical 
probability Pd at which a correlation mass diverges. From equation ([2]), (sq) diverges 
very slowly with respect to the system size as 

{so)l oc L ^ log2 N at p = p^. (3) 

Above Pel, {sq) diverges with the system size as 

{so)LOca^ = 2^^°^^" = N^. (4) 

where 

7/; = log2« = l-l/e. (5) 

The exponent ip mimics dj/d where df is the fractal dimension of a cluster embedded on 
the d-dimensional Euclidean lattice. Note that ip is a. function of ^ and increases with p. 
Even though a cluster mass diverges above Pd, the ratio {so)/N, which can be regarded 
as an order parameter, goes to zero for ^ cxd as far as < 1. li ^ diverges at a 
certain p, ip continuously approaches to 1. We define the critical probability Pc2 at which 
the correlation length diverges. At this critical point a prefactor of the connectedness 
function in equation ([1]) in form takes an important role. The exponent fj corresponds 
io 2 — d + rj for the d-dimensional Euclidean lattice. Positive fj makes {so)/N vanish in 
the thermodynamic limit and the order parameter continuously rising from zero at Pc2- 
In this case, however, {so)/N at Pc2 approaches to zero very slowly with A^ as (log2 N)'"^. 
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The finite size dependences are summarized as 

A^o for 

P < Pel 

log2 N at p = Pel 
(so) oc <j Ar'^(p) for Pci<p< Pc2 (6) 
iV/(log2iV)-^ at p = p,2 
N for p > Pc2- 

The cluster mass sliows fractal behavior in the phase Pd < p < Pc2 in contrast with the 
Euclidean lattice which has a critical point and unique fractal exponent. We consider 
the discrepancy of the probabilities for correlation mass and length is the key concept 
of the MPT on NAGs, which never occurs on the Euclidean lattices. In the following 
section we confirm the validity of the scenario shown above by numerical simulations. 

3. Numerical results 

Here let us explain the detail of our Monte-Carlo simulations. The systems with 
L = 10 — 22 are investigated. We generate samples, i.e., sets of open and closed 
bonds, by using pseudo random number generated by the Mersenne-Twister method 
[T6] . The ensemble average of observed quantities is taken over 480,000 samples. In 
order to improve the precision of the average values for small p, we treat the bonds 
which is connected to initial 3 generations exactly, i.e., counting all 2^^(^^~^)"^ possible 
realizations with probability — p)""", where Uo {ric) is the number of open (closed) 
bonds. This treatment only takes the CPU-time independent of A^. 

We also investigate the percolation on the dual lattice of the EBT (see figured]), 
which is also a NAG with one end. Each vertex of the dual lattice is put on the center 
of triangle or rectangle cells of the EBT and each bond, which is open with probability 
p, crosses with the conjugate bond of the EBT. A duality relation 

Pel + Pc2 = 1 and Pc2 + Pel = 1 (7) 

is known for dual planar NAGs |10j . 

As mentioned in the previous section the A^-dependence of (sq) for the finite size 
systems determines the expected three phases. So we estimate a fractal exponent 

rfln(so)7v d 

so that {sq) oc A^'^. In practice, we calculate ip by the difference [ln(so)2Ar — 
ln(so)Ar/2)]/[ln(2A^) — ln(A^/2)]. This is a good approximation of equation ([8]) for large 
A^ when (sq) is a power function of A^. The p-dependences of tp for both of the EBT 
and the dual lattice are shown in figure El Three distinct phases can be observed as 
expected. The exponent ip of the EBT grows from zero to unity in the intermediate 
phase between p ^ 0.30 and p ~ 0.56. Oppositely ip of the duality lattice decrease from 
unity to zero between 1 —p ^ 0.30 and 1 —p ^ 0.56. This suggests the duality relation 
actually holds. 
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Figure 2. The p-dependence of ip for the enhanced bmary tree and its dual lattice. For 
the latter, the horizontal axis is 1 — p to cofirm the duality relation. Results for various 
L's are shown together. The two vertical lines indicates p = 0.304 and p = 0.563, 
respectively. 



The A^- dependence of (sq) at Pd is decided as follows. Approaching pd from above, 
one expects that ip goes to zero but (sq) diverges as In TV (oc L) at p = pd- To exclude 
the possibihty that (sq) oc L*^ with ip 1, we plot 

as a function of p in figure [31 As L goes to infinity, v^l(p) decreases to zero for p < pd 
and diverges for p > pd- Only at p = pd, converges to A^-independent value, which is 
estimated at unity. Practically we obtain the precise critical probability Pd = 0.304(1) 
from this crossing point and Pd = 0.436(1) in the same way. Using the duality relation 
dB, Pc2 = 0.564(1) and p^^ = 0.696(1) are also determined. 

In figurelHwe show the geometry of the root cluster on the EBT by choosing samples 
which has a relatively large size cluster. Figure ID^a) is for p = 0.304 ^ pd- The cluster 
survive marginally; its branches do not show spreading behavior, and therefore the mass 
of the cluster is proportional to the number of generations. Figure ll](b) shows the cluster 
at a middle point in the intermediate phase, which has many branches. The branches 
are spreading with the generation and tend to diverge. However this cluster occupies 
very small part of the whole system because spreading rate is slower than that of the 
EBT itself. Consequently, there is a room for other infinitely large clusters. Figure 
mj^c) shows the cluster at p = 0.564 ^ pc2- It looks compact and will occupy the finite 
fraction of the whole system in the thermodynamic limit. 

Next we estimate the critical exponents of the first transition at p = pd- Let us 
suppose a finite size scaling law 

(So)ocL^^/%((Pel-p)i^^/'^), (10) 
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Figure 3. The p-dependence of ipL for the EBT (left) and its dual lattice (right) 
with several system sizes. The crossing point indicates the transition point where (sq) 
logarithmically diverges. 




for p < Pel- Q The scaling function so(-) should have asymptotic forms, 
„ , , f const, for X ^ 1 

So{x) (x\ (11) 
I X ' tor X ^ 1. 

FigureOshows a nice collapsing of data by using 7 = 1.0 and V/u = 1.0. The data for the 
dual lattice also shows same scaling behavior (the form of scaling function is very similar 
to the one for the EBT). This result confirms that the first transition at pd belongs to 
the mean field universality class (7 = 1 and V = v = 1/2 |T3| ) as predicted in Ref . [15] . 
On the other hand the transition at does not look like an usual continuous phase 
transition. The order parameter {sq)/N is fitted well by {so)/N = 0.49 + 0.58L~^'^^^ 
at p = 0.564 ^ Pc2 using data for L = 10 — 24 (not shown here). This means fj = 
{sq)/N has a finite limit value in L ^ 00 at the critical point. We have to be careful to 

X Consult Ref. [17| about the exponent 77. We use generation L as a characteristic length instead of 
chemical distance in these papers. The number of retrenched horizontal paths is so small compared 
with vertical ones that we regard this change to be irrelevant. 
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Figure 5. Finite size scaling plot of (sq) as a function oip — pci- The two groups, the 
EBT and the dual lattice, are plotted in arbitrary unit to collapse. 



conclude such special feature, discontinuous transition in the framework of second order 
transition, but can say, at least, that fj is smaller than 0.083. 

Finally we investigate the distribution function Ug, which is the number of clusters 
with size s per vertex, in order to show that the system is always critical in the 
intermediate phase as already implied from the A^-dependence of (sq). We assume a 
finite size scaling law 

n,{N) = N-'f'^nisN-'''), (12) 

for fixed p between p^i and Pc2- Here the scaling function n{x) is a power function x~'^ 
for X -C 1, and a rapidly decreasing function for x 1. In the scaling plot, figure El 
we use ip evaluated by equation ([8]) and the scaling relation r = 1 + ip~'^ at each p 
(note that this relation is based on the scaling relation r = 1 + df/d on the Euclidean 
lattices). Figure E] strongly supports that the above finite size scaling (IT^ holds in 
the intermediate phase. Equation f|T2l) indicates that the characteristic size of clusters 
diverges with N ^ oo and a power-law distribution, ns{s) oc s~'^, holds up to infinite 
s. When such a power-law exists, it can be said that the number of infinite clusters is 
infinite because one can always find tp' satisfying < ip' < ip for any positive tp with 
which the number 

POO 

N{s > N'^')= / dsNusiN) ~ iV(i-'/''/V')^ (13) 
diverges in the limit TV — oo. 

4. Summary and discussions 

In conclusion, we studied the bond percolation problem on the enhanced binary tree 
by Monte Carlo simulations for the first time. We observed the intermediate phase in 
addition to the nonpercolating and percolating phases where the average cluster size is 
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Figure 6. Finite size scaling plot of Us as a function of s for several p's in the 
intermediate phase. We use the values of ip shown in figure [21 Here we eliminate the 
data for (s < 16) since the data with too small s does not obey to the scaling law. 

finite and in the same order of tlie wliole system size A^, respectively. All the points in 
the intermediate phase is critical and the fractal exponent ip, where (sq) cx A^'^, increases 
continuously from zero to unity with increasing p. Note that the system looks critical 
in terms of mass of cluster, and the correlation length must diverge at Pc2- 

We expect that the mechanism of the MPT on the EBT is commonly observed in the 
other systems of NAGs with one end. There is, however, an interesting open problem, for 
both analytic and numerical studies, whether non-transitivity or inhomogeneity changes 
the scenario or not, which seriously corresponds to the complex networks. 
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